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Abstract 

We offer the parameter, interpreted as the "chemical potential", which is sen- 
sitive to the first order phase transition: it must decrease with number of evap- 
orating (produced) particles (hadrons) if the (interacting hadron or/and QCD 
plasma) medium is boiling and it increase if no phase transition occur. The main 
part of the paper is devoted to the question: how one can measure the "chemical 
potential" in the hadron inelastic processes. Our definition of this parameter is 
quite general but assumes that the hadron multiplicity is sufficiently large. The 
simple transparent phenomenological lattice gas model is considered for sake of 
clarity only. 
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1 Introduction 

^ . Despite the fact that the first order phase transition in the ion collisions is widely 

discussed both from theoretical [1] and experimental [2] points of view the feeling of 
some dissatisfaction nevertheless remain. To all appearance the main problem consists 
in absence of the single-meaning directly measurable ("order") parameter (s) which may 
confirm this phenomenon in the high energy experiment. Our aim is to offer such 
parameter, explain its physical meaning and to show how it can be measured. 

We guess that to observe first order phase transition it is necessary to consider 
very high multiplicity (VHM) processes. Then in this multiplicity region exist following 
parameter: 

(i,(n, s) ~- Tn,S ln<r w (s), (1.1) 
n 

Here a n is the normalized to unite multiplicity distribution which can be considered in 
the VHM region as the "partition function" of the equilibrium system, see Appendix, 
and T is the mean energy of produced particles, i.e. T is associated with temperature. 
Continuing the analogy with thermodynamics one can say that (— T In a n )/n is the Gibbs 
free energy per one particle. Then /i can be interpreted as the "chemical potential" 
measured with help of observed particles Q in a free state. 

We assume that the system obey the equilibrium condition, i.e. produced particles 
energy distribution can be described with high accuracy by Boltzmann exponent, or, it 
is the same, the inequality (b) must be satisfied. This assumption defines the "VHM 
region" [3] . It must be underlined that existence of the "good" parameter T does not 
assumes that the whole system is thermally equilibrium, i.e. the energy spectrum of 
unobserved particles may be arbitrary in our "inclusive" description. 

The definition (1.1) is quiet general. It can be used both for hadron-hadron and 
ion-ion collisions, both for low and high energies. It is model free and operates only 
with "external" directly measurable parameters. The single indispensable condition: 
we work in the VHM region of observed particles. It is evident that such generality has 
definite defect: measuring \x one can not say something about details of the process. 



1 Notice that one may consider n as the multiplicity in the experimentally observable range of phase 
space. 



This "defect" have following explanation. The point is that the classical theory of 
phase transitions have dealing immediately with the internal properties of media in 
which the transition occur. But in our, "S'-matrix", case one can examine only the 
external response on the phase transition in the form of created mass-shell particles. 

Continuing the analogy with the boiling, we are trying to define the boiling by the 
number of evaporating particles. The effect is evidently seen if the number of such 
particles is very large, i.e. in the VHM case. The "order parameter" is the work needed 
for one particle production, i.e. coincides with the "chemical potential". In the boiled 
"two-phase" region the media is unstable against "evaporation" of particles, i.e. the 
chemical potential must decrease with number of produced particles. 

We offer quantitative answers on the following three question. 

(A) In what case one may observe first order phase transition. 

We will argue that observation of VHM states are necessary to find phase transition 
phenomenon. First of all the energy of produced particles are small in VHM case. This 
means that the kinetic degrees of freedom does not play essential role, i.e. they can not 
destroy, wipe out, the phase transition phenomenon. The second reason is connected 
with observation that in the VHM region one may use such equilibrium thermodynamics 
parameters as the temperature T(n, s), chemical potential fi(n, s) and so on. 

(B) What we can measure. 

We will see that in VHM region exist the estimation (1.1) where all quantities in r.h.s. 
are measurables. 

(C) What kind effect one may expect. 

Chemical potential, /i(n, s), by definition is the work which is necessary to introduce, 
i.e. to produce, one particle into the system [I] . If the first order pase transition occurs 
then /i(n, s) must decrease with n in the two-phase ("boiling") region. It is our general 
conclusion which will be explained using lattice gas model. 



2 Definitions 
2.1 

We will start from simple generalization of well known formulae. Let us consider the 
generating function 

oo 

p(z, S ) = 5>V n (s). (2. 2) 

n=0 

For sake of simplicity a n is normalized so that 

p(l,s) = l. (2.3) 
One may use inverse Mellin transformation: 

dz 

2ni T z T 



ff »W = wra^' 4 ) (2 - 4) 



to find a n if p(z, s) is known. Noting that a n have sharp maximum over n near mean 
multiplicity n(s) one may calculate integral (mel) by saddle point method. The equation: 

d 

n = — — lnp(z,s) (2. 5) 

am z 

defines mostly essential value z = z(n, s). Notice that o~ n = if the hadron multiplicity 
n > n max = s/s/m, where m is the characteristic mass of hadron. Production of identical 



particles is considered for sake of simplicity. Therefore, only z < z max = z(n max , s) have 
the physical meaning. 

One may write p(z, s) in the form: 

p(z,s)=exp (2. 6) 

where the Mayer group coefficient hi can be expressed through fc-particle correlation 
function (binomial moments) Cfc(s): 

k=i K ' 



Let us assume now that in the sum: 

k\~ 



\np { z,s) = J2^-r ] —Ck(s) (2. 



k 

one may leave first term. Then it is easy to see that 

z(n, s) = n/ci(s), Ci(s)=n(s), (2. 9) 

are essential and in the VHM region 

ii 

InaJs) = -nln— -(1 + 0(1/ Inn)) = -n In z(n, s)(l + 0(l/lnra)). (2. 10) 
ci(s) 

Therefore, in considered case with Ck = 0, k > 1, exist following asymptotic estimation 
for n » 1: 

In a n ~ — n In z(ra, s), (2. 11) 

i.e. <7 n is defined in VHM region mainly by the solution of Eq.(eql) and the correction 
can not change this conclusion. It will be shown that this kind of estimation is hold for 
arbitrary asymptotics of a n . 

If we understand a n as the "partition function" in the VHM region then z is the 
activity usually introduced in statistical physics if the number of particles is not con- 
served. Correspondingly the chemical potential u is defined trough z: 

p = T\nz. (2. 12) 

Combining this definition with estimation (i.10) we define a n through «. But, if this 
estimation does not depend from the asymptotics of a n , it can be used for definition of 
fi(n, s) through cr n (s) and T(n, s). Just this idea is realized in (1.1). 



2.2 

Now we will make the important step. To put in a good order our intuition it is 
useful to consider p(z, s) as the nontrivial function of z. In statistical physics the 
thermodynamical limit is considered for this purpose. In our case the finiteness of 
energy and of the hadron mass m put obstacles on this way since the system of 
produced particles necessarily belongs to the energy-momentum surface^). But we can 

2 It must be noted that the canonical thermodynamic system belongs to the energy-momentum shell 
because of the energy exchange, i.e. interaction, with thermostat. The width of the shell is defined by 
the temperature. But in particle physics there is no thermostat and the physical system completely 
belongs to the energy momentum surface. 



continue theoretically a n to the range n > n max and consider p(z, s) as the nontrivial 
function of z. 

Let us consider the analog generating function which has the first n < n max coeffi- 
cient of expansion over z equal to a n and higher coefficients for n > n max are deduced 
from continuation of theoretical value of a n to n > n max . Then the inverse Mellin trans- 
formation (mel) gives a good estimation of a n through this generating function if the 
fluctuations near z(n, s) are Gaussian or, it is the same, if 

\2n — z 3 d 3 In p(z, s)/dz 3 \ 



\n + z 2 d 2 In p(z, s)/dz 2 \ 3 / 2 



« 1. (2. 13) 

z(n..s) 



Notice that if the estimation (1.10) is generally rightful then one can easily find that 
l.h.s. of (z) is ~ l/n 1 / 2 . Therefore, one may consider p(z,s) as the nontrivial function 
of z considering z(n, s) < z max if n max » n » 1. 

Then it is easily deduce that the asymptotics of a n (s) is defined by the leftmost 
singularity, z c , of in this way generalized function p(z, s) since, as it follows from (eql), 
the singularity "attracts" the solution z(n,s) in the VHM region. In result, we may 
classify asymptotics of a n in the VHM region if (z) is hold. 

Thus our problem is reduced to the definition of possible location of leftmost singu- 
larity of p(z, s) over z > 0. It must be stressed that the character of singularity is not 
important for definition of p(n,s) in the VHM region at least with 0(1/ Inn) accuracy. 
One may consider only three possibility: at n — > oo 

(A) z(n,s) -»• z c = 1, 

(B) z(n, s) — > z c , 1 < z c < oo, 

(C) z(n, s) — > z c — oo. 

Other possibilities are nonphysical or extremely rear. Correspondingly one may 
consider only three type of asymptotics in the VHM region: 

(A) o n > 0(e~ n ). We will see that in this case the isotropic momentum distribution 
must be observed, i.e. the energy, e, distribution in this case is Boltzmann-like, ~ e _/3e ; 

(B) a n = 0(e~ n ). Such asymptotics is typical for hard processes with large transverse 
momenta, like for jets [S] ; 

(C) o n < 0(e~ n ). This asymptotic behavior is typical for multiperipheral-like kine- 
matics, where the longitudinal momenta of produced particles are noticeably higher 
than the transverse ones [6] . 

We are forced to assume that the energy is sufficiently large, i.e. z max is sufficiently 
close to z c . In opposite case the singularity would not be "seen" on experiment. 

Our aim is to give physical interpretation of this three asymptotics. The idea, as it 
follows from previous discussion, is simple: one must explain the nature of singularity 
z c . It must be noted at the same time that in the equilibrium thermodynamics exist 
only two possibility, (A) and (C) [13] and just the case (A) corresponds to the first order 
phase transition. 

Summarizing the results we conclude: if the energy is sufficiently large, i.e. if z max 
is sufficiently close to z c , if the multiplicity is sufficiently large, so that (b) is satisfied 
and z(n, s) can be sufficiently close to z c , then one may have confident answer on the 
question: exist or not first order phase transition in hadron collisions. 

It must be noted here that the heavy ion collisions are the most candidates since z c 
is easer distinguishable in this case. 



2.3 



The temperature T is the another problem. The temperature is introduced usually using 
Kubo-Martin-Schwinger (KMS) periodic boundary conditions [7J . But this way assumes 



from the very beginning that the system (a) is equilibrium [8] and (b) is surrounded by 
thermostat through which the temperature is determined. The first condition (a) we 
take as the simplification which gives the equilibrium state where the time ordering in 
the particle production process is not important and therefore the time may be excluded 
from consideration. 

The second one (b) is the problem since there is no thermostat in particle physics. 
For this reason we introduce the temperature as the Lagrange multiplier (3 = 1/T of 
energy conservation law [3] . In such approach the condition that the system is in 
equilibrium with thermostat replaced by the condition that the fluctuations in vicinity 
of (3 are Gaussian. 

The interesting for us p(z, s) we define through inverse Laplace transform of p(z, (3): 

P{z,a)= ! jyLe*fip(z,P). (2. 14) 



It is known [8] that if the interaction radii is finite then the equation (of state): 

V~s = -^ P (z,P) (2.15) 

have real positive solution j3(n, s) at z — z(n, s). We will assume that the fluctuations 
near (3{n,s) are Gaussian. This means that the inequality [3] : 



\d 3 In p(z, (3) /d[3 3 



\d 2 \np(z,f3)/dp 2 \ 3 / 2 



« 1 (2. 16) 

z=z(n,s),(3=f3(n,s) 



is satisfied. Therefore, we prepare the formalism to find thermodynamic description of 
the processes of particle production assuming that this iS-matrix conditions of equilib- 
rium (z) and (b) are holcjj. 

We want to underline that our thermal equilibrium condition (b) have absolute 
meaning: if it is not satisfied then f3(n, s) loses every sense since the expansion in vicinity 
of /3(n, s) leads to the asymptotic series. In this case only the dynamical description of 
S'-matrix can be used. 

It is not hard to see [3] that 



* 1 



d 

■^InRlz, f3)\ z=z ^ s )^ = p M =< JJ(e~ < e >) >n, s (2. 17) 



i=l 



is the /-point energy correlator, where < ... > HiS means averaging over all events with 
given multiplicity and energy. Therefore (b) means "relaxation of Z-point correlations", 
/ > 2, measured in units of the dispersion of energy fluctuations, I = 2. One can note 
here the difference of our definition of thermal equilibrium from thermodynamical one 

m ■ 

2.4 

Let us consider now the estimation (1.1). It follows from (mel) that, up to the preex- 
ponential factor, 

lncr(n, s) ~ — n\nz(n, s) + ha. p(z(n, s), s). (2. 18) 



3 Introduction of P(n, s) allows to describe the system of large number of degrees of freedom in terms 
of single parameter, i.e. it is nothing but the useful trick. It is no way for this reason to identify entirely 
l/(3(n,s) with thermodynamic temperature where it has self-contained physical sense. 



We want to show that, in a vide range of n from VHM area, 



n In z(n, s) > In p(z(n, s), s). (2. 19) 

Let as consider now the mostly characteristic examples. 

(A) Singularity at z = 1. 
This case will be considered in Sec. 3. In the used lattice gas approximation In z(n) ~ rT 5 
and 

ln<7„ « n~ A = n\az{n){l + 0(l/n)). (2. 20) 

(B) Singularity at z = 1 + l/n 3 -(s) < oo: 

In s) = —7 ln(l — nj(s)(z — 1)). In this case z(n) = 1 + l/n^s) — l/nfij(s) and 

In cr n ~ —n/fhj(s) + 7 Inn = 

= -n/n i (s)(l + 0(lnn/n)). (2. 21) 

(C) Singularity at z = oo: \iap(z,s) = Ck(s)(z — l) fe , k > 1. 
In this case z(n) = (n/kck) 1 ^ >> 1 and 

lno- n « -n In z(n) (1 + 0(1/ Inn)). (2. 22) 

One can conclude: 

(i) The definition (1.1) in the VHM region is rightful since the correction falls down 
with n. On this stage we can give only the estimation of correction but (1.1) gives the 
correct n dependence. 

(ii) Activity z(n,s) tends to z c from the right in the case (A) and from the left if we 
have the case(B) or (C). 

(iii) The accuracy of estimation of the "chemical potential" (1.1) increase from (C) 
to (A). 



3 Ising model: phase transition 

The physical meaning of singularity over z [TD] may be illustrated by following simple 
model. As was mentioned above the singularity at z — 1 is interpreted as the first 
order phase transition. Therefore, let us assume [11] that (3 is so large that interacting 
particles strike together into clusters (drops). Then the Mayer's group coefficient for 
the cluster from I particle is 

6,09) ~ e ~^ l(d ~ 1)/d , 

where rl^W*, I >> 1, is the surface tension energy, d is the dimension. Therefore, if 
d > 1 the series over I in (2.4) diverges at al z > 1. At the same time, the sum (2.4) 
converge for z < 1. 

We consider following analog model to describe condensation phenomenon in the 
particle production processes. Let us cover the space around interaction point by the 
net assuming that if the particle hit the knot we have (—1) and (+1) in opposite case. 
This " lattice gas" model [4] has a good description in terms of Ising model [12] . We may 
regulate number of down oriented "spins", i.e. the number of produced particles, by 
external magnetic field Ji. Therefore the "activity" z = e _/3W , i.e. —Ti is the "chemical 
potential" [13] . 

Calculation of the partition function means summation over all spin configurations 
with constraint a 2 = 1. Here the ergodic hypothesis is used. It allows to exclude the 
time from consideration. 



To have the continuum model we may spread normally the ^-function of this con- 
straint 



5(a 2 - 1) ~ e 



-(1-<j 2 ) 2 /A 



Therefore, the grand partition function of the model in the continuum limit looks as 
follows Ej : 

p(ft,z) = [ Dae- S(a ' n \ (3. 23) 



where the action 



S(a,TC) = I d i x{-{Va) 



ea 2 + aa A — Act > . 



(3. 24) 



The structure of contributions in (2.6) essentially depends on the sign of constant e, see 
Fig.l where the case e > is shown. Following notations was used: 



(3. 25) 



where 1//3 C is the phase transition temperature. Phase transition takes place if (3 > /3 C 
(T < Tc), i.e. we will consider in present section e > 0. In this case the mean spin 
< ct >^ 0. We will assume that ft » ft c since in this case the fluctuation around < a > 
are small and calculations in this case became simpler. Considered model describes decay 
of unstable vacuum (181 . 




Figure 1: Solid line: undisturbed by TC potential and dotted line includes TC. 



The singularity over TC appears by following reason. At TC = the potential 

v = -ea 2 + aa A (3. 26) 

have two degenerate minima at ct = ±y^e/2a. The external field TC < we destroy 
this degeneracy. But in this case the system in the right-hand minimum (with the 
up-oriented spins) becomes unstable. 

The branch point in the complex plane corresponds to this instability. The discon- 
tinuity gives [H] : 

P(/M = ^#e-^ 2 , (3.27) 



where (/3 > f3 c ) 



vr2 (8pR 4 \ 7/2 f l-p c /p \ 3/4 R 4 

stt g / /? c y /2 i? 4 
j 81V2M p) imr 



(3. 28) 



It must be noted that the eqs. (eql) and (eq) have only one solution: 

(3(n, s) — > (3 C , In z(n, s) = l(n, s) — > 0. 

at increasing n. This means that the singularities at T = T c and z = 1 attracts the 
solution: 

l{n, 0) ~ n" 1 / 3 ^ - /3 C ) 7 / 6 , ^(n, s) = ^(1 + 7 /n 4 ), 



where 7 is the positive constant. 
In result, 

lnp n (s) ~ -n 2/3 (l/n 2/3 ) 7 ~ -1/n 4 ( nlnz(n, s)) 

decrees with n and the chemical potential 



■(1 + 7K) 



4\-l 



(3. 29) 



(3. 30) 



(3. 31) 



also decrees with n. 

Some comments will be useful to this Section: 

1. One may note that a n is defined by the discontinuity the the branch point in 
complex plane of In z and decay of the meta-stable states does not play any role. 




Figure 2: Stable ground state disturbed by Ti. 



2. It follows from (2.12) that, at fixed /3, 



\nz c ~ 



1/n) 1 / 3 « 1. 



(3. 32) 



This means that for large n our calculations are correct. At the same time, in VHM 
region z near unite is essential and it decrees with n. 

3. The work which is needed for production of one particle is ~ In z(n, s)//3(n, s). 
Therefore production of large number of particles needs less work per one particle. 



This conclusion have simple physical meaning (see beginning of present section). 
Let us consider decay of nonstable phase. The decay happens through production of 
clusters (domains with down oriented spins). The volume energy of cluster is oc i? 3 ,, 
where i? is radius of cluster. It burst the dimension of cluster. If R < R c , where R c 
is critical dimension of cluster, then the formation of such cluster is improbable. But 
if Rq > R c then the probability grows with radii of cluster. Latter explains why the 
chemical potential falls down with multiplicity. 

4. In the VHM region the temperature, T(n.s), tends to its critical value, T c , and 
slowly depends on n. 



4 Ising model: stable minimum 

Let as consider the system with stable vacuum, (3 < f3 c (T > T c ) in (2.8). In this case, 
see Fig.2, the potential v(a) has unique minimum at a — 0. Switching on external field 
H the minimum move and the average spin appears, (j(Ti) ^ 0. One can find it from 
the equation: 

—Act + 2ea + Aaa 3 = A, e > 0. 
Having d^Owe must expand the integral (2.6) near a: 

p(/3,z) = e Idx\a e -W{a)^ ^ 33) 

where W(a) expandable over a: 

00 p 1 

W ^ = J2l Uidx^n^B^,...^^ (4. 34) 

1=1 1 J i=i 

where Bi is the /-point one particle irreducible vertex function. In another wards, Bi 
play the role of virial coefficient. Comparing (3.2) with (2.4) one may consider a as the 
affective activity of group of / particles. 

The representation (3.2) is useful since in the VHM region the density of particles 
is large and the particles momentum is small. Then, remembering that the virial de- 
composition is equivalent of decomposition over specific volume, calculating Bi one may 
not go beyond the one-loop approximation, i.e. we may restricted by the semiclassical 
approximation. 

Therefore, having large density one may neglect the spacial fluctuations. In this case 
the integral (2.6) is reduced down to the the usual Cauchy integral: 

/ + 0O 
-oo 

In the VHM region A~7Y~ln,2>>lis essential. It is easy to see that 

a ~ -(A/4a) 1/3 (4. 36) 

is the extremum. The estimation of integral near this a looks as follows: 

p(/Moc{l2a(Ay /3 l ' e 3AVa /4(4Q)1 /3_ (437) 



4a 

This leads to increasing with n activity: 

l{n,s) ~n 8/3 (4. 38) 



and 

p(n,s) ~n 2/3 . (4. 39) 

In result, 

lnp n (s) n 11/3 (~ -nlnz(n,s)). (4. 40) 

A few comments at the end of this section: 

(i) Cross section falls dawn in considered case faster then 0(e~ n ). The estimation: 

Pn{s) ~ e - nln2(n ' s) 

gives the right expression in the VHM region. 

(ii) The chemical potential increase with n: 

^(n, s) = -(T(n, s) /n) In a n {s) ~ n 2 . (4. 41) 



5 Conclusions 

We may conclude that: 

(i) We found definition of chemical potential (1.1). This important observable can 
be measured on the experiment directly where T(n, s) is the mean energy of produced 
particles at given multiplicity n and energy ^/s. 

((ii) Being in the VHM region one may consider that fi(n, s) = 0(l/n) at compar- 
atively high multiplicities and it rise, fi(n,s) = 0(n) with rising multiplicity, Fig.3, at 
comparatively low multiplicities. The transition region is defines the critical tempera- 
ture T c . But it is quiet possible that the condition (b) allows to see only one branch of 
the curve shown on Fig.3. 

iv) The simplest example of finite z c presents the jet considered in the case (B), 
Fig.l. Hence case (C) has pure dynamical basis and can not be explained by equilibrium 
thermodynamics. 




T c n 

Figure 3: Chemical potential /i(n, s) as the function of multiplicity n. T c is the critical 
temperature. Breakthrough is the "two-phase" region. 
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